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Abstract. We show the preliminary results of our search for the progenitor systems of type la supernovae (SNe la). We model 
binary populations our aim being to compare these models with the observations of detailed element abundances of the hot 
Intra-Cluster Medium. 
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INTRODUCTION 

Type la supernovae (SNe la) are thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs). They are 
used as cosmological distance indicators [1, 2, 3, 4] and are the main source of iron in the Universe. However, the 
nature of their progenitors is not well understood. Two progenitor scenarios have been proposed, namely the single 
degenerate (SD) and double degenerate (DD) scenario. In the DD scenario, 2 CO WDs are formed close together by 
previous common envelope phases and merge due to gravitational radiation. In the SD scenario a CO WD accretes 
material from a companion that loses material by Roche lobe overflow (RLOF) or by a stellar wind. In this research we 
distinguish two types of companions, namely a hydrogen rich companion (SDn-scenario) and a more evolved helium 
rich companion (SDne-scenario). 

The observed delay time distribution (DTD) of SNe la, which is the rate of supernovae as a function of time since 
formation, can be used to constrain the progenitor models. In addition, the nucleosynthetic-yields of SNe la can be 
used to test the progenitor models. Clusters of galaxies are gravitationally bound systems in which metals from SNe, 
stellar winds, etc., accumulate over a long time span. They are ideal tools to examine the time integrated yields of SNe 
la. In this research we will study the evolution of SNe la using a binary population synthesis code based on single 
and binary evolution models [5, 6] in combination with a synthetic nucleosynthesis model [7, 8, 9]. We expect that the 
combination of the DTD, the chemical evolution of binaries and the observations of clusters will give us new insights 
into the evolution of stars towards type la SNe. 

BINARY POPULATION SYNTHESIS 

As a first step in this project we simulate binary populations in order to compare the resulting DTD to observations 
as well as to the results of other studies. The observed DTD shows a prompt peak [10] and a delayed component, 
which follows a power law, namely f^' [11]. We incorporate the efficiency of surface H and He burning on a WD 
and the model for an optically thick WD wind, based on Hachisu et al. [12] & Kato and Hachisu [13], into our 
binary population synthesis code. Other ingredients of our code and the initial distributions of binary parameters are 
as described by [5, 6]. In Fig. 1 we show the preliminary results of two simulations, in which we vary the common 
envelope (CE) efficiency; Oce = 3 (left panel) and oCce = 1 (right panel). The lower value of Oce results in closer 
binary orbits after the CE phase [6] . The simulations show that the three different progenitor channels contribute to the 
total rate of SNe la. Both our simulations show that the SDhc channel is dominant between 50 and 200 Myr, with an 
average rate over this time interval of 0. 1 -0.6 SNuM ' , depending on the common envelope prescription. The dominant 
scenario between 200 Myr and a Hubble time is the DD channel, with an average rate of 0.01-0.02 SNuM. The minor 
channel is SDh, with a rate about one magnitude lower then the DD channel in the case of a high Oce- The average 



SNuM = Supernova rate per 100 yr per lO'" M,;; in stars 




FIGURE 1. The Delay Time Distribution (DTD) of SNe la in SNuM of the three different progenitor scenarios in our model 
(Left Fig. acE = 3- Right Fig. oce = !)• The solid line indicates the rate of the SDhc channel, the dotted Une the SDh channel and 
the dashed line the DD channel. 



rate of type la SNe over a Hubble time is about 0.03 SNuM in both cases. Our results reproduce, at least qualitatively, 
the prompt peak and delayed dependence of the observed DTD. But our simulated rate is a factor of 2-3 lower than 
the observed one [11]. 

Other population synthesis studies also find a prompt peak caused by the SDhb channel [e.g. 14] and the power law 
dependence of the delayed SNe la [e.g. 15]. However, the largest deviation between the different models is in the SDh 
channel, due to the uncertainties in the accretion rate and accretion efficiency of H and He bunring on a WD and the 
efficiency of wind accretion onto a WD [12]. For example, [15] find 10~^ SNuM for the SDh channel, while [16] 
find a rate which is about an order of magnitude higher. 

With this research we hope to constrain the DTD by means of binary population synthesis and to ehminate the 
uncertainties in single and binary evolution. 
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